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Figure 1: Ligeia Mare (https://en.wikipedia.org/wiki/Ligeia_Mare) 


1 Introduction 


This application note simulates the speed of a small sea-going drone that is 
power surfing across Ligeia Mare (from Greek mythology). To start, it is 
assumed that the drone is moving with a uniform speed of three meters per 
second. After its electric engine is cut off, the drone begins coasting to a halt 
under the influence of atmospheric drag (assumed to be quadratic) and 
surface friction. The drone’s speed and the time it takes for it to come to a halt 
will be determined and displayed in figure 3. 


In section 2, a mathematical model will be developed for the drone’s speed as 
a function of time. 


For this application, an analytical solution, a numerical method solution, and 
an analog computation solution will be determined and then compared in 
table 1. 

2 Mathematical modeling 

Starting with Newton’s second law of motion, 


ma = Fpet = -drag - friction = -kv? — umg 


Assuming m = 1 kg, k= 0.1 kg/m, g = 1.352 m/s, up = 0.125 (just a guess), 
v(0) = 3 m/s, and neglecting units for simplicity, 


a= -0.1v? - 0.169 with v(0) =3 
Since a = dv/dt, 
dv/dt = -0.1v2 - 0.169 with v(0) = 3 (1) 


Equation (1) is separable and integrable. Additional mathematical details 
provided upon request: mcimorosi@desu.edu 


dv /dt = -0.1(v2 + 1.69) 


dv/dt = -0.1(v? + 1.32) 


Vv t 
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[eect fae 
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[tan?(v/1.3) - tan?(3/1.3)]/1.3 = -0.1t (2) 
v = 1.3tan[tan1(3/1.3) - 0.13t| (3) 


Setting v = 0 m/s in (2) to determine the time required to come to a halt, 


thait = tan-1(3/1.3)/0.13 = 8.94 seconds (4) 


3 Computer setup 
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Figure 2: Computer setup for coasting across the surface of a Ligeia Mare 
































Heun’s Numerical Method 
(hand-held programmable calculator) 


PROGRAM:TITANV 
:ClrHome:ClrDraw:Fix 2 
:"DV/DT=-0.1(V2 + 1.69)” 
“INITIAL CONDITION” 
(0-T:3-V 

:”STEP SIZE”:0.05-H 
:Lbl 1 

‘If V<-0.05:Then 

:Goto 2:Else 

:Disp {T,V} 

:-0.1(V2 +1.32)>F 
‘T+HOT 

:V+FH>W 

:-0.1(W2 +1.32)S 


(F+S)/2>A 


























:V+AH>V 
:Pause:Gotol 
:Lb]12:End 
4 Results 
t(s) Analog Analytical Solution | Heun’s Numerical Method 
Computation v(m/s) Solution H = 0.05 s 
Solution v(m/s) 
v(m/s) 
0.00 3.00 3.00 3.00 
2.50 1.50 1.44 1.44 
5.00 0.80 0.73 0.73 
7.50 0.30 0.25 0.25 
9.50 0.00 
8.94 0.00 
8.95 0.00 




















Table 1: Solution Comparisons 





Figure 3: Drone speed as a function of time* 


*For this application note, the display was produced during a single run by a 
differential equation analog computer prototype using discrete components 
with tolerances between 1% and 10%. 
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